arXiv: 1505.07464v2 [hep-ph] 25 Sep 2015 


Dark matter ignition of type la supernovae 

Joseph Bramante^ 

^Department of Physics, University of Notre Dame, 225 Nieuwland Hall, Notre Dame, IN, USA 

Recent studies of low redshift type la supernovae (SNIa) indicate that half explode from less 
than Chandrasekhar mass white dwarfs, implying ignition must proceed from something besides 
the canonical criticality of Chandrasekhar mass SNIa progenitors. We show that 1 — 100 PeV mass 
asymmetric dark matter, with imminently detectable nucleon scattering interactions, can accumulate 
to the point of self-gravitation in a white dwarf and collapse, shedding gravitational potential energy 
by scattering off nuclei, thereby heating the white dwarf and igniting the flame front that precedes 
SNIa. We combine data on SNIa masses with data on the ages of SNIa-adjacent stars. This 
combination reveals a 2.8cr inverse correlation between SNIa masses and ignition ages, which could 
result from increased capture of dark matter in 1.4 versus 1.1 solar mass white dwarfs. Future studies 
of SNIa in galactic centers will provide additional tests of dark-matter-induced type la ignition. 
Remarkably, both bosonic and fermionic SNIa-igniting dark matter also resolve the missing pulsar 
problem by forming black holes in > 10 Myr old pulsars at the center of the Milky Way. 


It would be difficult to overstate the significance of type 
la supernovae, which have been used to measure the vac¬ 
uum energy of the universe. The uniform width-to-height 
ratio of SNIa emission curves, known as the Phillips re¬ 
lation, allows SNIa distance and redshift to be inferred 
from luminosity over ^40 day outbursts [T]. A preva¬ 
lent lore states that the uniformity of SNIa light curve 
width-to-height ratios results from a uniform population 
of SNIa progenitors: carbon-oxygen white dwarfs (WDs) 
explode after reaching Chandraskehar mass {Mch 
1.4 Mq) by accretion from a binary companion. 

However, a recent multiband study of 20 “la-norm” su¬ 
pernovae light curves, selected from 147 low redshift spec¬ 
imens, presents strong evidence that many type la super¬ 
novae do not reach Mch before exploding [2]. Follow-on 
studies have concluded that about half of the observed 
SNIa have sub-Mcn progenitors [3], resulting in a rela¬ 
tion for SNIa progenitor masses, 

MsNia/MQ = 1.322 ± 0.022 -b (0.185 ± 0.018)a:i, (1) 


ence for heavier (1.4 Mq) WDs exploding sooner, which 
is one prediction of SNIa-igniting DM. 

In the analysis of DM-induced WD ignition that fol¬ 
lows, it is sufficient to assume any model with some 
amount of DM asymmetry (we consider totally asym¬ 
metric DM, but partly asymmetric DM only changes 
the effective capture rate) and a mass around a PeV. 
These requirements allow enough DM to accumulate, 
self-gravitate, and collapse within a WD lifetime. We 
also assume a velocity-independent DM-nucleon cross- 
section, (Tnx- Note that, because the DM momentum 
transfer to standard model (SM) particles is less than a 
GeV throughout, this simplified framework can be UV 
completed with the addition of weak scale mediators. 
PeV mass particles can fill out the DM relic abundance 
(despite unitarity bounds (TU]) through, for example, 
non-thermal processes m or CO- annihilation enhanced 
freeze-out [T3] . 

1. DM accumulation. DM’s collection rate in a WD 



where xi parameterizes SNIa light curve stretch and 
ranges from —2.8 to 2.6. Clearly, these findings are at 
odds with Mc//-attainment as the sole mechanism for 
triggering SNIa. 

There are proposed mechanisms for triggering sub- 
Chandrasekhar SNIa with known particles, e.g. accreted 
helium shells can cause “double-detonations” ms] or bi¬ 
nary WD mergers [7]. These require binary companions, 
and lone WD SNIa progenitors are preferred by some ob¬ 
servations, including the absence of expected luminosity 
shocks from binary companions [5], little circumstellar 
material in progenitor systems (corroborated by low ini¬ 
tial x-ray and radio emission), and the high overall rate 
of SNIa [3] . 

In this paper we demonstrate that heavy asymmetric 
dark matter (DM) ignites lone WDs by rapidly increasing 
the WD temperature inside the region of DM collapse. 
This occurs for halo density DM {px ~ GeV/cm^) col¬ 
lected into a WD within its lifetime. We correlate SNIa 
masses with nearby stellar ages and find a 2.8cr prefer¬ 


( 2 ) 

where G is Newton’s constant, px is the DM halo den¬ 
sity, Myj and are the WD mass and radius, mx is 
the DM mass, v ^ 200 km/s is the WD-DM velocity 
dispersion, a sat ~ Rw/^n is the maximum DM-nuclear 
cross-section {Nx is the number of nuclei in the WD). We 
take h = c = ks = i throughout. The square-bracketed 
term accounts for DM that scatters but is not captured 
in the WD, and = QmxVssc/i^Nv'^i'm'X/’mx — 1)^, 
where Vesc — \/‘^GMys /Rw is the escape velocity from 
the WD surface and mAr ~ 14 GeV is the average WD 
nuclear mass. The next sections discuss the coherent 
scattering form factor incorporated in the DM-nucleus 
cross-section, cjax- (If the DM-WD momentum trans¬ 
fer exceeds the inverse nuclear radius (1/3 fm“^) during 
capture, we calculate DM capture on nucleons instead of 
nuclei.) 
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Minimum fi„x prompting DM collapse and SNIa ignition 
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Figure 1. The top panel displays the minimum DM-nucleon 
cross-section necessary for a WD to accumulate a collapsing 
sphere of DM in t™ = 0.5 Gyr (thin cyan lines) and tm = 5 
Gyr (thick red lines). The galactic DM density assumed to 
surround WDs is px ~ GeV/cm®. Solid, dashed, dotted- 
dashed and dotted lines correspond to WD masses, radii, and 
core WD densities indicated. The bottom panel shows the 
minimum DM-nucleon cross-section required for asymmetric 
DM to ignite SNIa after collecting and collapsing in WD pro¬ 
genitors with central temperatures of Tui = 10^ K. For lighter 
DM, the number of DM particles required for collapse can be 
too large - the collection rate is limited by the WDs geo¬ 
metric cross-section - leading to the low mass cutoffs shown. 
For heavier DM, the number of DM particles required for col¬ 
lapse decreases, along with the rate of energy injected into the 
WD, leading to a high mass cutoff for WD ignition. These 
calculations assume that the DM-DM scattering cross-section 
equals the DM-nucleon cross-section axx ~ f^nx, that is 
we assume the mediator responsible for DM-nucleon scatter¬ 
ing induces comparable DM-DM scattering. Each panel also 
indicates that heavy asymmetric fermionic (or bosonic with 
quartic A ~ 1) DM collapses <Gyr old pulsars in the Milky 
Way’s galactic center, assuming a central parsec DM density 
px^ ~ 10'* GeV/cm®. The LUX spin-independent 2a bound 
on DM-nucleon scattering is given in purple |18[ llOj . 
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After the DM is captured, it will continue to scatter 
with the WD until it thermalizes at its center. DM with 
a mass lighter than 100 PeV and cross section anx ^ 
10“*^ cm^ will thermalize in less than a hundred million 
years [20] . The thermalized DM sphere will have a radius 


given by the Virial theorem, 

rth = {^TwI^T^Gpw'mxY^'^ 



where = Pu,/(10® g/cm®), and = lO'^- 10® g/cm® 
is the density of a sub-Chandrasekhar, ignitable WD. If 
the DM particles are heavier, the thermal radius will be 
smaller, and fewer DM particles need collect within a 
WD lifetime to become self-gravitating and collapse. The 
amount of DM required for collapse is 

Nsg = Axp^r^f^/imx ^2x 10 ®® 



In Figure we show the DM-nucleon cross-section re¬ 
quired for DM to collapse inside = 0.9 — 1.4 Mq 
WDs within 0.5 and 5 gigayears. 

2. DM collapse. As the sphere of DM collapses, it 
will shed gravitational potential energy by scattering off 
carbon and oxygen nuclei in the WD, which can prompt 
SNIa. The rate and mode of DM collapse is determined 
by whichever of the gravitational free-fall time and the 
DM-DM interaction time is shorter m- At the onset 
of collapse in a WD, weakly-interacting DM will not be 
self-thermalized. As collapse progresses, the dynamical 
free-fall time (which scales with the DM sphere’s radius 
as oc r®/®) will come to exceed the DM-DM scattering 
time (oc Once the free-fall time exceeds the DM- 

DM interaction time, the DM will self-thermalize through 
DM-DM interactions, and (in the absence of a more effi¬ 
cient radiation mechanism) shed gravitational energy by 
nuclear scattering as it collapses further. We will see that 
self-thermalized, collapsing DM heats the WD enough to 
spark SNIa. 

In total, there are three relevant timescales for the 
collapsing DM sphere: the dynamical free-fall time, the 
DM-DM interaction time, and the DM-nuclear interac¬ 
tion time. The free-fall time is 


iff 


V37r/(32Gpx) 0.15 s 



( 5 ) 


Note that at the time of collapse, the DM density equals 
the WD density. 

The DM-DM interaction time is 


txx 


(nxfTxxvx) ^ ~ 3 X 10® s (Pws) ^ 


Tw \ / o'xx 

10^ Ky V10“*® cm^ 



( 6 ) 


where nx is the DM number density and vx = 
^j2GNsgmx/r is the DM velocity for a self-gravitating 
DM sphere of radius r. 
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The DM-nucleus interaction time is 

taX ~ {nN<yaXVx)~^ — 10^ S (Pws)”^ 

where the final expression includes the DM-nuclear form 
factor (at maximum coherence), which we discuss shortly. 
Comparing the dependence on the collapsing sphere’s ra¬ 
dius (r) in Eqs. ([^, and Q, we see that heavy DM 
in WDs begins collapse in free-fall, but will transition to 
thermalized collapse once the DM-DM interaction time 
becomes shorter than tfj. Setting tff = txx, we solve 
for the radius below which the collapsing DM is thermal¬ 
ized by DM-DM interactions, 


^sta 


V NsgCJxx - 0.1 cm 



107 Ky V10-40 cm2 / 


( 8 ) 


Prior to collapsing to Vsta^ the DM sphere must discard 
the difference in gravitational potential energy between 
rth and Vsta by scattering off WD nuclei. Setting r = rth 
in Eq. Q, and noting that a DM particle will reach 
equilibrium with the WD after ~ mx/^mx scatters pT] . 
we estimate the maximum time required for the DM to 
collapse to radius Tgta, 


/ mx \ 0/2 

tth < iaxmx/^iJiM ^ 0.1 yrs (Pws) 

( ( anx \-4 

^VlO^xy V10-40 cm2 / 


( 9 ) 


In the analysis that follows, we assume (t„x ~ <^xx, that 
is, that the same mediator responsible for DM-nucleon 
scattering also mediates DM-DM scattering, giving both 
processes about the same cross-section. Relaxing this 
assumption, particularly axx ^ cTnX, results in other 
DM-igniting SNIa parameter space, as we will discuss. 

3. SNIa ignition. Once the DM has self-thermalizes 
(at Tsta), it uniformly heats the WD (within rgta) as 
it collapses further. First note that the DM-WD mo¬ 
mentum transfer is smaller than the inverse of the WD 
nuclear radii, p ~ mxVsta < (3 fm)-^ where Vsta = 
yj2GNsgmx/rsta- (It can also be verihed that Vsta is the 
relative velocity of the DM-WD nuclear system.) This 
indicates that scattering off WD nuclei will be coher¬ 
ently enhanced. The DM-nucleus cross-section is Gax — 
A^(3ji[a;]/a;)^Exp(—x^/3)cr„x, where = 200 accounts 
for coherent scattering enhancement off carbon and oxy¬ 
gen nuclei, ji is the Bessel function of the hrst kind, 
X = pvn, and r„ ~ 3 fm [211 122 ] ■ The average energy 
transferred per DM-nuclei scatter is then e ~ mxVg^^j2, 
occurring on a timescale set by tsta = taxirsta)- The 


rate of energy transferred to the WD is 


Qhe = Nsge/tsta 2 X 10^^ GeV/s 

\107 K J V10-40 cm2/ V10-40 cm2 


- 23/8 





where this assumes coherent DM-nuclei scatters, but we 
use the preceding nuclear form factors in computations. 

If Qhe is larger than the rate at which the heat dif¬ 
fuses in the WD, the collapsing DM can ignite SNIa. 
The work of [25 showed that for WD material of mass 
irihe, in the mass range 10’ ® < (mhe/g) < 104®, SNIa 
ignition requires heating a mass mhe of carbon-oxygen 
to temperature (T/^e/lO® ’^ K)70/3 > (pu,8)^^^(g/wfte)- A 
close inspection of Eq. Q reveals that heavier DM in 
cooler WDs will enclose less than a milligram of WD 
at Tsta- Therefore, to remain well within the numeri¬ 
cal calculations of [23|, we consider heat diffusion out 
of a sphere of radius r/jg, where is the larger of 
Tsta and a sphere enclosing a milligram of WD, = 
1.5 X 10-4 cm (pws)”^^^- Then for SNIa ignition, the 
DM heat transfer (Eq. @) must exceed the conductive 
diffusion rate out of a r/jg size sphere at the WD center, 
given by [22] 


Qdif An'^rheTheiThe - T^)/15 kcPw, (H) 


where Kg ~ (IQ-^ cm^/g)(The/107 K)^'®(pw8)~^'® is the 
conductive opacity for WDs, which have a thermal diffu¬ 
sion dominated by relativistic electron conduction when 
Pw ^ 10® g/cm® [ 22 ] . The preceding expression for Kg 
conforms to WD conductive opacity tables over the rele¬ 
vant range of WD temperatures and densities [2312S] ■ 

In Figure [l] we display DM masses and minimum 
cross-sections required to ignite WDs within 0.5 and 5 
Gyr. In addition to Qhe > Qdif, we require that the 
DM impart enough heat to bring the WD core tem¬ 
perature to The- This requirement is easily fulfilled 
if Qhe > Qdif, because the capacitance of a degener¬ 
ate WD’s ion lattice is simple (c« = 3fc;,/2 [22])) so in 
the limiting case of slowest heating {mx ~ 100 PeV, 
p-ui — 10® g/cm®, and cr„x — 10-42 cm®), the DM 
will heat a gram of WD to lO^® K (1 MeV) after a 
fraction of collapsing DM particles have scattered once, 
Nsca ~ (10®® nuclei/g)(MeV/e/Vsg) < IQ-^. The curves 
in Fig. [^cut off at high mass, because too few DM par¬ 
ticles collapse, c.f. Eq. 0. to adequately heat the WD. 

It can be shown that PeV mass fermionic DM is not 
degenerate (and is not Pauli-blocked) as it collapses 
through Tsta- To calculate the radius at which the DM 
becomes degenerate, first note that the most energetic 
stabilized fermions have minimum Fermi kinetic energy 
Ef = (^xNsgjjTmxT^■ This implies a DM radius 
(using the Virial theorem). 


Tdeg = (97r/4)®/®/Gm|Vy/® - lO’^ cm 


X 


107 K 


- 1/2 



{PWif^^, 


( 12 ) 
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SNIa progenitor age versus mass 

(MsNla/Mo) 



Figure 2. Data bins (with la error bars) are given for SNIa 
age versus progenitor mass, by combining the SNIa-adjacent 
star age data in [5S] with the progenitor mass fitting func¬ 
tion of [3]. Asymmetric DM model curves for SNIa-ignition 
age as a function of progenitor mass are overlaid (these were 
obtained using the same methods as Figure [^. Note that 
the DM model curves assume uniformly carbon-oxygen WDs 
with temperature 10^ K. 

smaller than by two orders of magnitude. A similar 
computation reveals that collapsing, heavy bosonic DM 
will not condense until r <C rsta- 

On the other hand, if axx ^ o^nX (we note again 
that this analysis makes the simplifying assumption 
that axx ^ cTnx), the radius at which the DM self- 
thermalizes will be smaller. In this limit, we have found 
that fermionic DM collapsing into a degenerate sphere 
(and bosonic DM condensing into a Bose-Einstein con¬ 
densate m), will heat and ignite SNIa for much of the 
parameter space indicated in Figure 

4. SNIa age vs. mass. Results in prior sections 
indicate that if SNIa are triggered by PeV mass DM, this 
implies a negative correlation between SNIa progenitor 
age and mass. In Figure we reframe a study that 
correlated host galaxy star age with light curve stretch 
pS] (see also [HHSI]), by converting light curve stretch 
to SNIa mass with results from 0 (see Eq. [^. Refs. [2H] 
and [3] use SNIa stretch parameters s and xi respectively; 
Refs. [5^ l551 convert s to Xi. 

In Figure^ 67 SNIa data points from Ref. [55] are 
collected in bins spanning 0.5 — 5 and 5—14 Gyr. Our 
analysis of vertical error bar heights agrees with [28] . 
The vertical separation between the two mass bins, ac¬ 
counting for uncertainty introduced converting s to Xi 
to MsNia, is 0.15 M© {Mbini = 1.42 ± 0.035 Mq, 
Mbin 2 = 1.27 ± 0.041 Mq), amounting to a 2.8a sig¬ 
nificant separation, comparable to the 3.3cr result in |28j . 
which correlates age with s instead of SNIa mass. We 
overlay model curves for SNIa-igniting asymmetric DM. 

5. Galactic center pulsar implosions. The cen¬ 
ter of the Milky Way does not harbor as many pulsars 
as expected [M] [3S] . Asymmetric DM, more dense in the 
galactic center (GC), could abundantly collect in GC pul¬ 
sars and form pulsar-destroying black holes, as first noted 
in [3^. DM models fitting the missing pulsar anomaly 


were presented in [ 37 ] [ 35 ], and this putative population 
of imploding pulsars could also be the source of fast radio 
bursts [ 33 ]. Figure shows that PeV mass asymmetric 
DM would destroy 0.1 Gyr old pulsars at the GG, where 
we use the same pulsar calculations as [33 [ 35 ]. 

The pulsar-imploding parameter space in Figure [l] is 
insensitive to anx , because the pulsar DM capture cross- 
section saturates when anX ~ 10“"^® cm^, meaning the 
amount of DM collected in GC pulsars is constant for 
o'nX ^ 10“"^® cm^. With the mass of collected DM re¬ 
maining constant, heavier DM fermions will form black 
holes in pulsars while lighter DM will not, because the 
critical mass necessary for black hole formation drops 
as the DM mass increases, ~ Mpjm\. PeV 

mass bosonic DM with an order one quartic self-coupling, 
A ~ I, will have the same critical mass for forming a black 
hole as fermionic DM, ~ [JD]. This is 

particularly important when comparing these results to 
studies that assume asymmetric bosonic DM with a van¬ 
ishing quartic (A < 10“^®) [27l l4lll47] . Note also that 
the cross-sections we consider are smaller than those rel¬ 
evant for other stellar probes of DM (e.g. main-sequence 
[48H53] . astroseismic [54H56] . cooling [57H6T] . and pulsar 
phase [62] constraints). 

6. Conclusions. We have introduced a mechanism 
for igniting SNIa. Heavy asymmetric DM with de¬ 
tectable SM interactions can collapse in and heat 0.9-1.4 
solar mass WDs, prompting a thermonuclear runaway. 
We have pointed out a possible inverse correlation 
between SNIa progenitor masses and ages, which is 
predicted if DM ignites SNIa. Studies of SNIa in galactic 
centers, which harbor a denser bath of DM and would 
prompt younger and less massive WDs to explode, 
will provide an additional probe of DM-ignited SNIa. 

Note added. As this paper was being finalized. Ref. 
[53] appeared; it also considers ways, mostly different, 
that dark matter can trigger supernovae. 
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